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Abstract 

£h A 3D-2D dimension reduction for a nonlinear optimal design problem with a perimeter penalization is performed in 
i_ijhe realm of T-convergence, providing an integral representation for the limit functional. 

I Resume 

^ Reduction dimensionnelle 3D— 2D d'un probleme non lineaire d'optimisation de forme avec penalisation 
sur le perimetre On effectue dans ce travail une reduction dimensionnelle 3D-2D d'un probleme non lineaire d'op- 
.timisation de forme avec une penalisation du perimetre. Une representation integrate de la fonctionnelle limite est 
obtenue. 

y—\ ' 

(N- 

Version francaise abregee 

x- ,. 

5^ On s'interesse dans ce travail au comportement asymptotique d'une suite de problemes non lineaires d'optimisation de 
5^ forme avec penalisation du perimetre sur le domaine cylindrique 51 (e) := ui x (— s, e) , ou e > et ui est un ouvert borne 
de M 2 . On suppose que le domaine fl (e), occupe par le solide, est constitue par deux materiaux hyperelastiques dont les 
densites d'energie, Wi : K 3x3 — > R, i = 1, 2, sont continues, verifiant la condition de croissance 

& (\F\ P - 1) < Wi(F) < /3(1 + \F\ P ) pour tout F G R 3x3 , p>\, i = 1,2, avec > > 0. (1) 

Plus precisement, on considere le probleme de minimisation suivant 



inf {-( f (XB{e)W l + {l- X E{e))W2){Vv)dx- f f e ■ V <Lr 

XE {e) e W(Q(e);{0,l}) (2) 

+P(E(e);Q,(e))} : v = sur du x (-e,e), C 3^ £ ^ J Xe{ s ) dx = a|, 

f2(e) 

ou E (e) C fi (e) est un ensemble mesurable de perimetre fini et / £ G L p (r2(e);R 3 ) oil ^ + ^ = 1 et A G [0, 1] est la 
fraction du volume rempli par le premier materiau. 

Nous commencons par effectuer un changement de variables afin de rendre le domaine independant de e (cf. (7)). On 
obtient ainsi le probleme suivant 



inf 

ueW^ p (n;R 3 ) 
X G BV(SI;{0,1}) 



{ X W 1 + (1 - X )W 2 ) ( V qU |-V 3 u )dx- I / ■ ucfe 



(ft): 



(3) 



= sur da; x (-1, 1), ^j^y J Xdx = a| . 



En utilisant les techniques de T— convergence, on demontre que, lorsquc e tend vers zero, le probleme (3) converge vers 
un probleme nonlincairc bidimcnsionnel. 

II suit un resultat general ou la penalisation du perimetre dans lc probleme initial est remplacee par une integrale 
elliptique. Ainsi, on considere \P : IR. 3 — > [0, +oo[ une fonction paire, continue, positivement homogene de degre 1 et telle 
que 



3C e]0, +oo[ : W e R 3 —\u\ < < C\v\. 
On etudie done, le comportement asymptotique, lorsque e tend vers zero, du probleme suivant 

inf {!( f {xE(e)W 1 + (l-XE(e))W 2 )(Vv)dx- [ f £ -vdx 

»eW^(fi( E );l 3 ) £ a { e 



(4) 



(<0 



fi(e) 



X6B7(fi( £ );{fl,l}) 



+ y ^(^(e))^ 2 ) : v = sur 9^ x (-e,e), £3^^ ^ Xe( £ ) tfo = a|, 

3B(s) fl(e) 



(5) 



ou "H 2 designe la mesure de Hausdorff restreinte a dE(e) et Ve( s ) es t la normale exterieure a E(e). On obtient finalement, 
le probleme limite (19). 



1. Introduction and setting of the problem 

The study of thin structures has been the object of many investigations. In particular, in mechanical engineering it is 
important for applications to minimize, under a given system of loads, the compliance (namely, the opposite of the total 
energy at equilibrium) of a given structure, satisfying a constraint on the volume. In order to design thin structures 
with the best possible resistance- weight ratio, the asymptotic behaviour of the compliance as the thickness of the sample 
tends to zero is studied. For a background on the modelling of thin plates we refer to the monographs of [1] and [11]. 

Let (e) := lo x (—£,£) , where to is a bounded open domain of M 2 and s > 0, and for the sake of illustration let 
us assume that Q(e) is clamped on its lateral boundary. We suppose also that 0(e) is hlled with two materials with 
respective energy densities W\ and W2, where Wi : R 3x3 — >• R, i = 1,2 are continuous functions satisfying (1) (see 
Remark 2 below, where assumption (1) is discussed). 

Let E(e) be the first phase, f e be the given load on ft(e) and assume that the volume fraction of each phase is given by 
^ : ~ £3(Q( £ )) Jn( £ ) XE(s){ x )dx G [0, 1], where Xe(e) denotes the characteristic function of the phase E{e). The compliance 
C £ (xE(e)) 1S defined as 

C% XE{S) ) := - yewi M (£);R3) {1 J ((xe^W, + (1 - XE{S) )W 2 ) (Vv) -f s .v)dx:v = on x (s, s)}. 

fi(e) 

In [6] , [7] and [8] the asymptotic behaviour of a 3D optimal elastic compliance problem is studied, as the thickness (or 
the cross section in the case of beams) tends to zero and the volume fraction in the design region remains unchanged. 
It is assumed that the material has a convex and 2-homogeneous potential and the analysis is performed in the small- 
displacement setting. The prescription of the volume in the minimum problem can be dropped by adding a Lagrange 
multiplier to penalize the volume in the cost functional. The asymptotic analysis performed in these papers leads to 
a fictitious material with local density, taking all the values in [0, 1] and not to a precise limit set, due to the loss of 
compactness in the characteristic functions. 

In this work we focus our attention on studying the worst possible design of a two-phase mixture of elastic materials 
in a thin film in the same spirit of [13] and [9], where the asymptotic analysis of a two-field minimization problem has 
been studied (i.e. (x, u), (design region, deformation)) as the thickness of the sample tends to zero. Having in mind the 
results contained in [2] and [17], we introduce a perimeter penalization in our functional in order to derive from the 3D 



Let us consider the optimal design problem in (2) where E (e) C £1 (e) is a measurable subset of 0(e) with finite 
perimeter, i.e., 



P{E(e);Sl(e)) := sup j ^ diupda: : tp G C*(fi(e);R 3 ), |M| L ~> < l} < +00, 



(6) 



E(e) 



and we assume that the load f £ G -L p '(Jl(e); M 3 ), with A + ^ = 1. 

In order to study the asymptotic behavior of (2) we reformulate our problem in a fixed 3D domain through a -- 
dilation in the transverse direction x 3 and then we perform T— convergence with respect to the pair (design region, 
deformation). Set 



:= oj x (-1,1), F £ := {(xi,x 2 ,x 3 ) G : (xi,x 2 ,ex 3 ) G £(e)}, 

u (xi,x 2 ,x 3 ) := w (xi,X2,ex 3 ) , / (xi,x 2 ,x 3 ) := f £ (xi,x 2 ,ex 3 ) , xe e {xi,x 2 ,x 3 ) := \E(e) (xi,x 2 ,ex 3 ) , 



(7) 



where v is any admissible field for (2). 

In the sequel we will denote x a := (x\,x 2 ), dx a := dx\dx 2 and V a and D a will be identified with the pair 
(Vi,V 2 ), (Di,D 2 ), respectively. For every matrix F G M. 3x2 and any z G K 3 , F := (F\z) represents the matrix 



p3x3 



whose first two columns are those of F and the last column is given by the vector z. 



Observe that by (6) and using the definition of total variation, P (E (e) ; Q (e)) = \Dxe(e) \ ( e ))- Making the change 
of variables y 3 := ex 3 and y a := x a we have A |-Dx.E(e)| (^( £ )) = | (D a Xe \^D 3 Xe)\ (^)> where Xe : = X£ e stands for the 
characteristic function of F e . Hence we are led to the rescaled minimum problem (3). 

For every e > 0, let J £ : L 1 ^; {0, 1}) x L p (Vt; M 3 ) -> [0, +00] be the functional defined as follows 



Je(X,u) := < 



XW1 y a u 



V 3 « + (1 - X)W 2 V a u 



-V 3 u^J^Jdx — J f-udx- 









\d 3X ) 







(SI) 



in BV(Sl; {0, 1}) x W hp (n; K 3 ), 
otherwise. 



(8) 



+00 

Let V : {0, 1} x R 3x3 -> [0, +00) be given by 
:= X W 1 (F) + (l-x)W 2 (F), 
with Wi and W 2 satisfying (1). Analogously, let V : {0, 1} x M 3x2 -> [0, +00) be defined as 

V* (x, ^) : = XW 7 ! (^) + (1 - X) "^2 (F) , with W l (F) := irrf Wi (F|c) , F G M 3x2 , i = 1,2. 



(9) 
(10) 



Consider the functional J : L^fi; {0, 1}) x L p (tt; 



[0, +00] as 



Jo(x, u ) : = ' 



1 

2 / QV( X , V a u)dx a - J Jf- udx a dx 3 + 2\D aX \(u), if (x,«) G BV(w; {0, 1}) x W^'^R 3 ), 



■i» (11) 
+00 otherwise, 
where QV stands for the quasiconvexification of V in the second variable. Namely, for every (%, F) G {0, 1} x R 3x2 



QV (x, F) := inf j J V (x, F + V aV ) dx a : ip G (Q'- » 3 ^ 



(12) 



where Q' C R 2 denotes the unit cube. 
Our main result is the following. 

Theorem 1.1 The family of junctionals { J e } I '-converges, with respect to the strong topology o/L 1 (fi; {0, 1}) xL p (Sl; } 
to Jq, as e — > + . 

Remark 1 We observe that Theorem 1.1 entails the convergence, as e — > + , of problems (2) in their rescaled version 
(3) to the problem 



inf \ 2 J QV(x, V a u)dx a - j j f ' ■ udx a dx 3 + 2\D a x\{w) : ^ J X dx a = \ X \- 



(13) 



In fact, this is due to the strong convergence in Z/ 1 (fl; {0, 1}) x L' p (Vl: R 3 ), of the sequence of almost minimizers {(xs, u £ )} 
of (3) to (x,u) G BV(u>; {0, 1}) x W ' p (u>;M. 3 ). And so, the constraint in the volume fraction C i^ E ^ Jn( £ ) XE(e)dx = 
£ j(Q) JsiXedx — A is kept in the limit, as well as the boundary conditions (cf. Remark 2). 

It is worthwhile to compare Theorem 1.1 with similar results in [13] and [9]. To this aim we observe that, in spite 
of what is proven therein, namely T-convergence results with respect to the convergence £^L ak>t x L p for (%, u), the 
presence of the perimeter in our energy (8), allows us to have a stronger convergence on the characteristic functions and 
thus to determine the worst possible design set. On the other hand, the fact that the perimeter is inserted in our model 
leads naturally to compare our results with those contained in [2]. Indeed, if W\ and Wi are of type W\ (•) := a' |-| 2 and 
W2 (•) := a \-\ , with < a 1 < a suitable constants, clearly QV(x, F) coincides with a'x|F| 2 + a(l — x)\F\ 2 - Hence, by 
[2, Theorem 2.2] the solution of the minimum problem (13) is locally Holder continuous and the optimal design set is 
equivalent to an open set A x (—1, 1), A C w. More refined results, about regularity in 2D, in the convex setting, may 
be found in [14,15] and in the references quoted therein. 

In our model, the lack of convexity in W\ and Wi entails, as underlined by Proposition 2.2, that we obtain a limit 
energy which depends continuously on the characteristic function of the design set and requires a quasiconvexification 
procedure in the deformation variable. We refer to [10] for regularity results related to our setting. 

Details about the results are contained in the next section, while for the properties related to T- convergence, sets of 
finite perimeter and BV functions we refer to [12] and [3], respectively. 



2. The limit problem 



We start by stating the properties of the energy densities in (9) and (10) that we will exploit in the sequel. 
Proposition 2.1 Let V be as in (10) . Then V is continuous and satisfies 



P'(jF\ P -l)<V(x,F)<f3(l + \F\ P ), 



where j3' and j3 are the constants in (1) . Moreover, \V (%, F) - V(x', F)\ < 2/J |x - x'\ (1 + \F\ P ). 
Proposition 2.2 The function QV in (12) is continuous and satisfies (14), an d 

\QV( X ,F) - QV( X ',F)\ < C\ X ' - x|(l + \F\ P )- 



(14) 



(15) 



Remark 2 We claim that energy bounded sequences {(xe,u e )} for problem (3), with u e clamped on duj x (—1,1), are 
compact m L^O^O, 1}) x L p (ft;R 3 ) and with limit in L 1 (w; {0, 1}) x L p (uj;R 3 ). 

If {{Xe, u e)} is a sequence such that J e (xe,u e ) < C, then there exists C G M + such that the following bounds hold 











-V 3 u e 




(DaXe 






LP 





(fi) < c. 



An argument entirely similar to that exploited in [16, Lemma 3], entails that there exists u G H^'P^R 3 ) such that 
V3U = 0, and so u can be identified with a function (still denoted in the same way) u G W 1,p (to; R 3 ). Thus we may find 
a subsequence, not relabelled, {u £ } such that u e — u in W 1,p (ft; R 3 ), and a measurable set E C fl such that Xe^ * Xe 
and D$xe = 0. Hence, there exists E 1 c 10, with \Dxe\(Q) — 2\Dxe> |(w), where E — E' x (—1, 1). In the following we 
will identify the set E with the set E' and denote xe' by x- 

Proof of Theorem 1.1 For every e > 0, let J e be the functional in (8). The T-convergence with respect to the separable 
metric space L 1 ^; {0, 1}) x L p (f2;R 3 ) ensures that for each sequence {e} there exists a subsequence, still denoted by 
{£}, such that T - lim^o+^^KO, 1}) x L p (n;R 3 ))J e exists. 

For every (x, u) G L 1 ^; {0, 1}) x L p (ft; R 3 ), let J (x, u) be its T-limit. By virtue of Urysohn property, it suffices to prove 
that any sequence {J e } admits a further subsequence whose r-limit, J(x, u), coincides with J (x,u) in (11). 

We observe that if (x, u) G (L 1 (f2; {0, 1}) x L p {Vt- R 3 )) \ {BV(u; {0, 1}) x W^ p (lu; R 3 )), then J(x, u) = +00. Indeed, 
if this is not the case, from J(x, u) < +00 we would get the existence of a sequence {(xe, u e )} converging to (x, u) such 
that J E (xe,u E ) < +00 and by Remark 2 this would imply (x,u) G BV(u; {0, 1}) x W 1,p (lj; R 3 ). 

The remaining proof is divided into two steps. First we show the lower bound, then we prove the upper bound. 
Step one: We claim that for every (x, u) G BV(u; {0, 1}) x VF^ p (o;;] 



J(x,u) >2 / QV (x (x a ) , V a u (x a )) dx a 



f (x a ,x 3 ) u (x a ) dx a dx 3 + 2 \D a x\ (w) . 

To prove the claim, let {(x e ,M e )} C L 1 (J7;{0, 1}) x L P (£!;R 3 ) be a sequence converging to (x, u ) S BV(ui; {0,1}) x 
W 1,p (w;R 3 ). For the forces and for the perimeter the lower bound follows by LP strong convergence of {u e } and lower 
semicontinuity of the perimeter, respectively. 



scaled gradients {(V a w £ , ^V 3 w £ )} are p-equiintegrable, and C 3 (fl\A £ ) — > as e — > + , where A £ := {x G il : u £ = w £ }. 
Denoting the bulk energy density of J £ by V as in (9), one obtains 



liminf / V [ Xe, ( V Q u e 



dx 



> liminf / V ( \e, ( V Q w £ 



> liminf 

£^0 + 



Xe, ^ a W £ 



-V 3 w £ 

e 

-V 3 w £ 

e 



dx — (3 lim sup 



1 



V a w £ 



-V 3 w £ 

e 



dx 



(16) 



n\A e 

dx > liminf / V (Xs, V Q w e ) dx > liminf / QV (% e , V a w e ) dx. 

e^0+ J £ ^0+ 7 

n n 



Observe that, by (15) 

f \QV(xe,V a w e )~QV( X ,V a w £ )\dx<C J \xs-x\(l + N a w £ ndx. 



(17) 



Thus, the p-equiintegrability of {(V Q w £ | jrV3W e )} and (17) ensure that as e —¥ + , Xe can be replaced by x m the 
right-hand side of (16). 

The density QV(x{x a ), •) is quasiconvex in M 3x2 for a.e. x E ft. Using an argument similar to that exploited in [16, 
Proposition 6] one concludes that QV(x(x a ), •) is quasiconvex also in M 3x3 . Thus, by the growth condition of QV, 
as stated in Proposition 2.2, the functional v G W 1,p (f2;R 3 ) i — > J n QV(x{x a ),V a v(x))dx is sequentially weakly lower 
semicontinuous with respect to H /1 ' p -weak topology. Hence, 

liminf / QV(x, V a w £ )dx > 2 / QV(x, V a u)dx a . 



By the superadditivity of the lim inf we achieve the claim. 

Step two: To prove the reverse inequality we start by observing that, fixing \ G BV(uj; {0, 1}), 
J{x,u) < liminf e ^ + J £ {x,u £ ) for every {u £ } C LP(f};R 3 ) and u G W^wjR 3 ) such that u £ -> u in L*>(ft; 
Thus it suffices to study the asymptotic behaviour with respect to the M^ 1,p -weak convergence of 



xWi V a u 



V 3 u £ + (1 - x)W 2 V Q u, 



i V 3 u £ ^dx- J f ■ u £ dx. 



Since x 1& fixed we can rewrite xWi(-) + (1 — X)W2{-) as a new function with explicit dependence on x a . 
Namely, let W : u x R 3x3 -> R be given by 

W(x a ,F) := V( X (x a ), F) = xMW^F) + (1 - x (x a )) W 2 (F), 

for every (x a , F) G cj x R 3x3 . 

Clearly, 14 7 is a Caratheodory function satisfying the growth condition ^j\F\ p — C < W(x a , F) < C(l + \F\ P ) for a.e. x a G 
w and for all F G R 3x3 . Applying [4, Theorem 2.3] to the sequence of functionals {G £ }, where G £ : L p (Vl\ R 3 ) -> [0, +oo) 

1. 



is given by G £ (u) := < 



n 

+oo 



V 3 m rfa; 



/■uda; if u G W 1 ' p (£l; R 3 ), 
otherwise, 



and arguing as in [9, Remark 3.3] we get 



lim 

E-S-0+ 



W [x a ,[ V Q M, 



\ V 3 u £ ))dx-jf- u £ dx\ < 2 / V a u)dx a - / / / • udx a dx 3 , 

n / oj -1 u 



where W : uj x R 3x2 is defined by W(x a ,F) := inf cgR 3 W(x a , (F\c)), and stands for the quasiconvexification of 14 7 
in the second variable. 
Observing that by (10) 

W(x a ,F) = x{x a )W[{F) + (1 - x{x a ))W 2 {F) = V( X (x a ),F), and QW{x a ,F) = QV( X (x a ),F) 
for every (x Q ,.F) G oj x R 3x2 , the proof is concluded. 

In the following we apply the previous analysis, with small changes, to the case where the perimeter penalization in (2) 



Recall problem (5) and observe that by J dE ^ ^(i>E(s))drl 2 we mean the integral with respect to the Hausdorff measure 
concentrated on dE(e), identified with S XB .., with exterior normal Ve(s)- 

By performing the same rescaling as in (7) we obtain the following formulation of (5) in the fixed domain SI, 



inf 

u e W 1>p (n;R 3 ) 
X e BV(SI;{0,1}) 



( X W 1 + (1- X )W 2 ) V Q 



-V 3 u 

e 



dx 



f ■ u dx 



+ 



J * (^l^s) dV. 2 : u = on duj x (-1, 1) , -^^y J X dx = a}, 



(18) 



where X denotes the characteristic function of E e , and v the normal to its jump set, S x . 

Next we will identify functions defined in R 3 (or R 2 ) and their restrictions to S 2 (or S 1 ), and so the same notations 
will be adopted. 

Let * : R 2 -> [0, +oo) be the function given by ¥(77) := inf {#(77, £) : £ G M} , with * : R 3 ->■ [0, +oo) as in (4). 
Consider the following minimum problem 



inf 

X &BV(u;{0,l}) 



QV(x,V a u)dx c 



f ■ udx a dx 3 + 2 1 * (u^dH 1 : 



1 



-1 u> 



Xdx a 



1 



> . 



(19) 



where ^P** denotes the convex envelope of Namely, *** (v) := sup{.g : R 2 — > R : 5 is convex < ^(w) Vu € R 2 }. 
Theorem 2.1 Consider the problems (18), and their minimizers. Then the latter converge, with respect to the strong 
topology o/L 1 (il; {0, 1}) x L p (il;M. 3 ), to the minimum of problem (19). 

We conclude by observing that the density "J" satisfies all the well established properties for the lower semicontinuity 
of surface integrals, such as BF-ellipticity (cf. [3, Definition 5.13 and Theorem 5.14]), since any continuous even function 
4> : S N ~ X — > [0, +00) is By-elliptic if and only if its positive 1-homogeneous extension is convex. 
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